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I. INTRODUCTION 


In this letter we wish to confirm and discuss in detail new data on 
the total omnidirectional cosmic y-ray flux in the 0.25 MeV to 6 MeV range 
obtained by the Environmental Research Satellite- 18 (ERS-18) which were 
presented in a preliminary form earlier (Vette, et al. , 1969). Diffuse 
cosmic y-rays between 0. 07 MeV and 1 MeV were first detected on the 
Ranger III (Metzger, et al. , 1964) and confirmed between 20 keV and 
0. 5 MeV by balloon and satellite measureme nts (Bleeker and Deerenberg, 
1969; Schwartz, 1969). Extension of the spectrum to higher energies has 
become of considerable importance since the recent successful detection 
of galactic y-rays >100 MeV by the OSO-III (Clark, Garmire and Kraushaar, 
1968). The ERS-18 measurements indicate an additional component of 
y-rays in the 1-6 MeV range above that expected by simple extrapolation 
of the spectrum below 1 MeV to higher energies. 

II. DESCRIPTION OF EXPERIMENT 
The ERS-18 was launched April 28, 1967, along with the Vela 4 
satellites into a highly elliptical orbit of 117, 500 km apogee and 15, 000 km 
perigee. The 7. 8 Kg satellite (Peterson et al. , 1968a) contains a number 
of radiation detectors designed to measure magnetospheric particles, 
cosmic-rays, solar X-rays and cosmic y-rays. A large detector, designed 
to measure the cosmic y-ray spectrum in the energy range 0.25 - 6 MeV, 
consisted of 7. 65 cm long x 6. 35 cm diameter Nal(Tl) scintillation counter 
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with a 1 cm thick plastic anti- coincidence shield for charged particle 
rejection. The shield counter covered all but one circular face of the Nal 
crystal. Six channels of pulse height analysis permitted spectral measure- 
ments in the five energy loss bands given in Table I. In addition, total 
counting rates were monitored in integral channels at 25 keV and 6 MeV. 

The anti- coincidence was switched on and off by the spacecraft commutator 
about every 500 seconds in order to assess the counting rates due to cosmic- 
ray energy losses. The difference in counting rate for these two modes 
served as a calibration using the cosmic- ray beam and verified operation of 
the anti- coincidence system. 

The satellite spent a large fraction of its 48-hour orbit in regions 

where magnetospheric particles were not detectable by the charged particle 

detectors whose lowest threshold energies were 40 keV for electrons and 

380 keV for protons. Data reported here were obtained primarily during 

May and June, 1967, when the local time of satellite apogee was generally 

directed toward the dusk meridian. Near apogee at 18 earth radii the earth 

-2 

subtends a solid angle of about 10 steradian; earth albedo y-rays therefore 
produced a completely negligible flux. Furthermore, because of the high 
perigee (2.4 earth radii) the satellite did not appreciably penetrate the high 
energy trapped proton zones which could induce radioactivity in the satellite 
and the Nal crystal. This effect has plagued previous experiments on the 
OSO-I (Peterson, 1965) and on the ERS-17 (Peterson, 1968b) but did not 
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contribute an observable background rate to this experiment. Furthermore, 
the data presented here were selected from intervals during which there was 
no evidence of detectable particle fluxes in interplanetary space above those 
attributable to galactic cosmic rays. Data obtained during the solar proton 
events of late May and early June have been excluded from the analysis. 
Although the information bandwidth of the ERS-18 telemetry was 6 Hz, and 
only a small fraction of this was devoted to the y-ray experiment, the 
measurement was nevertheless limited by systematic effects rather than 
counting statistics. 


III. DATA ANALYSIS 

Rates obtained with the anti- coincidence on and off are shown in 

Table I and Fig. 1, after converting to energy loss spectra by dividing by 

2 

the channel width and by the isotropic geometry factor, 54 cm . Since the 
difference between the two rates must be due to cosmic-ray effects, the 
data shows that energy losses due to cosmic- rays and photons are about 
equal in the highest channels. The rates are not a true measure of the 
incident spectrum since neither photons nor cosmic- rays necessarily lose 
all their energy upon interacting with the detector. Obtaining the actual 
input spectrum is a difficult process and generally requires some knowledge 
of the photon spectrum to energies well above the highest energy loss channel. 

The expected energy loss spectrum of minimum ionizing cosmic- rays 
isotropic on the detector, including edge effects, has been computed using a 


4 


Monte Carlo method and an interplanetary cosmic- ray intensity of 
2 

3. 07 /cm -sec which was measured on the Vela satellites simultaneously 

(J. R. Asbridge, private communication). The computed spectrum, shown 

as a dotted line in Fig. 1, is flat at energy losses below that corresponding 

to ionization losses by particles traversing the diameter of the crystal. 

That this calculated spectrum due to particle losses agrees so closely with 

the observed one is taken as substantial evidence that the anti- coincidence 

performed properly and that the energy loss calibration had not appreciably 

changed in orbit. Furthermore, the total integral rate (>6 MeV) measured 

in the central detector, after a 10 per cent edge effect correction, agrees 

with the independently measured cosmic-ray flux by Asbridge to about 5 per 

cent. This provides additional confidence in the entire system operation. 

The energy loss spectrum with the anti- coincidence on must 

therefore be due to y-rays incident on the detector from space or produced 

in the satellite itself. No known process has sufficient cross-section for 

cosmic- rays to produce an appreciable y-ray flux in passing through the 
2 

8 gm/cm of average satellite material adjacent to the y-ray counter. 
Although a detailed Monte Carlo calculation of the cascade development in 
a small amount of matter has never been carried out, one can estimate 
that only about 5 per cent of the cosmic rays will undergo nuclear inter- 
action in traversing the satellite. The electromagnetic and nucleonic 
cascades which follow develop to their maximum intensity after about 
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100 gm/cm . Our direct experience on balloon flights indicated that 

. 2 

7 gm/cm of A1 surrounding detectors similar to the ERS-18 does not 
result in an additional flux detectable above cosmic- ray produced atmo- 
spheric y- rays. 

A further experimental measure of the coupling between cosmic-rays 
and measured y-rays can be obtained by searching for correlations of the 
Y-ray counting rates with the modulation of galactic cosmic- rays by solar 
activity associated with the events of May- June 1967. During this period 
the sea level cosmic-ray indices varied some 10 per cent and the intensity 
at the satellite by 30 per cent. The correlation coefficients between the 
integral (> 6 MeV) channel and the 1-2, 2-3, 7, and 3. 7 - 6. 0 MeV channels 
were typically 0. 3 with the anti-coincidence on and 0. 5 with the anti- 
coincidence off. We have also averaged the counting rate in each y-ray 
channel over two specific time periods--one when the > 6 MeV channel was 
counting less than 142 c/sec and one when it was counting this value or 
greater. The differences between these averages were small compared 
with the differences between the averages for the > 6 MeV channel or the 
neutron monitor indices over the same time period. From these results 
we can infer that the contribution from y-rays produced in the satellite is 
less than 10 percent +_ 10 per cent for the highest energy channels. The 
comparison of our cosmic-ray rates with the Vela results, as well as 
estimates of background effects, lead us to assign an estimated total 
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systematic error to these fluxes of not more than 20 per cent. The relative 
error between the channel-by- channel results are considerable smaller 
than this. 

IV INTERPRETATION OF RESULTS 
The ERS-18 data are shown in Fig. 2 compared with previous 
results obtained on Ranger III in space (Metzger, et al. , 1964) and on the 
OSO-I in low altitude earth orbit (Peterson, 1966). These measurements 
were all obtained with detection systems similar to the ERS-18, and none 
provided direct information on the directionality of the incident radiation 
except for the inference that the incident y-rays cannot be highly anisotropic. 

A portion of the Ranger III counting rate above 1 MeV attributed to an 
instrument artifact is not shown. The OSO-I upper limit is obtained directly 
from the counting rate over the geomagnetic equator using the 2. 8tt solid 
angle factor subtended by space at 550 Km as the conversion factor. The 
interpretation of the OSO-I data was considerably more complicated than 
those of the ERS-18 and Ranger III because of background effects of (a) the 
2 00 Kg spacecraft, (b) earth albedo y-rays, (c) induced radioactivity and 
(d) gain calibration problems. These effects produce a measured total flux 
considerably above the ERS-18 result. However, from the distinct latitude 
dependence of the counting rates, it seems clear that an appreciable fraction 
of the OSO-I rate must be due to cosmic-ray production. An OSO-I lower 
limit based on the assumption that the background latitude effect due to earth 
albedo and spacecraft production is less than that due to primary cosmic-rays 


is also shown in Fig. 2. 
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Since the difference in the measured cosmic- ray flux between the 
ERS-18 and Vela 4 is essentially zero, a restriction on the integral y-ray 
spectrum may be made assuming that the maximum cosmic y-ray contribu- 
tion must be less than that due to the estimated systematic errors in the 

ERS data. The twenty per cent maximum errr r corresponds to a > 6 MeV 

. 2 

photon-produced rate of 0. toO c/cm -sec; the true integral rate is perhaps 
much less. 

Indicated in Fig. 3 and Table I are theoretical fluxes based on the 
interpretation of the 0. 01-1.0 MeV y-rays as due to Compton scattering 
of intergalactic electrons on the 3° K radiation. These spectra are normal- 
ized in the 30 - 300 keV range where good agreement is obtained among the 
various experiments. The original idea of Felten and Morrison (1966) 
indicated a simple power law spectrum whose differential number index 
was about -2. 3. The more recent spectrum of Brecher and Morrison (1969) 
takes into account the distribution of radio galaxy spectral indices and 
luminosities as well as cosmological factors. At energies above a few 
hundred keV, the energy loss spectrum as measured by the detector cannot 
be directly compared with either of the two theoretical photon spectra. The 
theoretical spectra which are fitted to the data below 1 MeV after correction 
to energy loss will produce counting rate spectra well below the measured 
data points in the 1-6 MeV range. We interpret this as indicating an 
additional cosmic y-ray component above that accounted for by extrapolating 
the known diffuse background to higher energies. No reasonable extrapolation 
of lower energy data can completely account for the ERS-18 data. 
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In order to find the spectrum of the excess component, Monte Carlo 

computer programs have been developed to compute the energy loss 

spectrum in the detector for any incident photon spectrum. The calculations 

for the two theoretical diffuse background spectra shown in Fig. 3 may be 

compared with the measured spectrum and differences taken to find the 

excess rates. The results appear in Table I. Data in the 0. 25 - 0. 6 MeV 

range have not been included in the following analysis because of (a) the uncertain 

correction for satellite absorption, (b) the rather wide energy window, and 

(c) the near agreement of the two theoretical spectra at lower energies. 

The total excess counting rate above that accountable by the extrapolated 

2 

diffuse component varies from 0. 31 to 0. 36 c/cm -sec over the 0. 6 - 6. 0 
MeV range and is 2. 5 to 5. 5 times larger than the extrapolated diffuse back- 
ground, depending on the model spectrum used to extrapolate from lower 
energies into this range. 


V. DISCUSSION 

These results must be discussed in terms of other measurements 
of the total cosmic Y-ray spectrum. Based on an analysis of balloon data, 
Bleeker and Deerenberg (1969) have found the spectrum over the 20 - 220 keV 
range to be even steeper than we have used with an index 2.45 + 0. 10. The 
cosmic flux in the 7. 9 - 110 keV range is known to be isotropic to within a 
few per cent (Schwartz, 1969) except for the contributions due to galactic 
point sources. No directional measurements are available in the 0. 2 - 10 
MeV range. The OSO-UI >100 MeV results indicate both a diffuse 
background component and a galactic "line" source. 
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We assume that the ERS-18 y-ray flux contains a component due 

to an extrapolation of the isotropic diffuse background spectrum and an 

2 

excess counting rate of about 0. 35 c/cm -sec in the 0. 6 - 6 MeV range. 

Since nothing is known about the directionality of the excess component, 
several possible origins must be considered. 

a) Galactic Point Sources . A collection of point X-ray- type 

sources radiating strongly in the MeV range could produce the excess. 

Although no X-rays greater than 500 keV from known sources have yet been 

-3 2 

detected, firm upper limits on the Crab Nebula of 1. 7 x 10 c/cm -sec-MeV 

_3 2 

at 1 MeV and 1. 0 x 10 c/cm -sec-MeV at 5 MeV have been obtained 
(Peterson, et al. , 1966; Schwartz, private communication). These limits 
are factors of 1. 3 and 2. 5 above an extrapolation of the hard X-ray spectrum 
of the Crab Nebula to higher energies. To account for the excess flux 
observed by ERS-18 would require some 50 sources radiating at the present 
Crab upper limit and even more using its extrapolated X-ray spectrum. At 
least some of such a large number of discrete sources radiating this strongly 
would most likely have been discovered by now. 

b) A Galactic "line" Component. A galactic disc component of 

-4 2 

cosmic y-ray radiation with an average intensity 2. 0 x 10 /cm -sec- rad 

has been detected at 100 MeV (Clark, Garmire and Kraushaar, 1968). 

2 

Explanation of our excess of 0. 31 c/cm -sec as a line component would 

2 

require an average flux of . 020 c/cm -sec-rad in the 1-2 MeV range. 
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Extending the high energy spectrum steeply to lower energies is ruled out 
by the work of Clark, Garmire and Kraushaar (1968) who state that the 
line source spectrum is rather flat in the 100 MeV range. In order to be 
consistent with the upper limits obtained on the OSO-III over the 7. 7 - 210 keV 
range (Schwartz, 1969), explanation of the ERS-18 excess in terms of a 
galactic "line" source requires a new component which is flat below 1 MeV 
and falls off steeply above 6 MeV. 


c) An Extragalactic Component . If the ERS-18 excess is to be 
interpreted as an additional isotropic component above that due to the 
spectrum extrapolated from lower energies, it seems clear this component 
must fall steeply in the 6-100 MeV range to connect with OSO-III 100 MeV 
results. The OSO-III experiment indicates, in fact, that the diffuse cosmic 
source spectrum at 100 MeV is steeper than that of the galactic line source. 

A number of explanations have already been advanced to explain the 
ERS-18 results. The spectrum of Brecher and Morrison (1969), which fits 
the data in the 3 keV to 1 MeV range rather well, also apparently fits the 
measured ERS-18 energy-loss spectrum. However, as already shown in 
Table I, when proper account is taken of the detector response this source 
fails considerably to account for the new results.* 

Clayton and his co-workers (Clayton, Colgate and Fishman, 1969; 


Clayton and Silk, 1969) have suggested that y-ray fluxes in the 1-3 MeV 

, j ....56 6. ld_ 56 77d 56 56 

range may be produced by the process Ni — »Co — »Fe if Ni is 


It also conflicts with the OSO-III datum, which was plotted a factor of ten too 
high in their paper. 
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synthesized in supernovae explosions. Each decay chain produces an 
average of nearly 5 y-rays, with energies extending to 3.47 MeV. Clayton 
and Silk have integrated the fluxes expected due to a reasonable rate of 
supernovae explosions per galaxy over cosmological distances and predict 
a detectable flux at the earth. This explanation requires a rapid fall off in 
the flux above 3.26 MeV; this was not observed on the ERS-18. Furthermore, 
they also predict 77d decay times to be occasionally observed at the higher 
energies due to supernovae in nearby galaxies. Although the total variation 
observed during the May- June period was only 5 per cent, the nearly 
continuous data extending over a 13 -month lifetime of the ERS-18 is only 
now being reduced to look for 77-day effects. Although this particular 
nuclear process does not explain the new data, the possibility of nuclear 
y-rays cannot be excluded. 

Stecker (1969a, 1969b,) has proposed that tt° meson decay 
photons occurring as a result of cosmic-ray interactions at a high density 
epoch of the expanding universe would be red shifted and result in 
a spectrum peaked in the 1-10 MeV range for a ^ max > the earliest time 
of cosmic-ray acceleration, of about 100. Since Stecker’s theory makes 
definite spectral predictions, we have attempted to fit his models to the new 
data in the 0. 6 - 6. 0 MeV range making maximum use of the ERS-18 energy 
loss spectrum. By using the normalization as a minimization parameter, the 
energy loss spectrum computed from the Stecker photon spectrum of a given 
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Z is least- square fitted (the relative errors are minimized) to the two 
max 

excess energy loss spectra given in Table I and shown again for convenience 

in Table II. The results of this fitting are indicated in Table II and Fig. 3 

using a Z of 50, 70, 100 and 150. The total photon spectrum lies a 

max 

factor of 1. 5 to 3, depending on the energy, above the energy-loss spectrum 
and clearly indicates the essential importance of correcting for the detector 
response. The standard deviation of the relative errors between the fit 
and measured values, a re j , are also given in Table II. The errors 
assigned to the calculated values are based on least square fitting differ- 
ences rather than the smaller relative errors in th.e measured data. 

A comparison of calculated and observed excess indicates that 

Z =100 provides the best fit for a power law background, and Z =70 
max max 

results in the best fit for the Brecher-Morrison spectra. The photon fluxes 

for these cases are shown plotted in Fig. 3. Although the spectra based on 

the Stecker model can be adjusted to fit the ERS-18 excess in the 0. 6 - 6. 0 

range rather well, the predicted flux at 100 MeV lies at least a factor of three 

- 6 2 

above the measurement of 1.2x10 photons /cm - sec-Mev- ster (letting 
A E = E) reported from the OSO-III measurement at 100 MeV. Of course, 
if one uses the data point at 100 MeV in conjunction with the ERS-18 data, 
a higher Z value will be obtained which produces a poorer fit for the 
0. 6 - 6. 0 MeV interval. Since the reported 100 MeV results may be high 
by a factor of two or three (Garmire, private communication), it is prema- 
ture to pr^rent such a fit. Furthermore, the prediction of the Stecker model 
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at 100 MeV depends mainly on the assumed cosmic ray spectra’ ndex at Z^x 
rather than on Z saax itself. A determination of 2 ma y is best made near the 
peak of the gamma ray spectrum. 

We believe the ERS-18 data have indicated clearly an additional 
component of the cosmic Y-ray spectrum. The Stecker model provides a 
photon spectrum which quantitatively agrees with the measured results by 
using a Z max somewhere between 70 - 100. Several other possible source 
mechanisms cannot be definitely ruled out by the present experiment. New 
experiments with greater spectral resolution, higher energy coverage, and 
directional sensitivity are needed to test theoretical models more definitively. 
The important cosmological implications of these models should give considerable 
impetus to carrying out such experiments. 
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photons /cm - sec 


TABLE II 


CALCULATED AND OBSERVED ERS-18 RATES USING STECKER MODEL 




Power Law Background Subtracted 
Rate (counts /cm^- sec- MeV) 


Energy Loss 
Range (MeV) 

Z = 50 

a . = .26 
rel 

Z = 70 

• 13 

Z = 100 

a . = . 096 
rel 

Z = 150 

a = . 2 1 

rel 

Obs. 

Excess 

0. 6 - 1.0 

0. 125 + . 032 

0. 166 + . 022 

0. 202 + .019 

0. 244 + . 051 

0. 191 

o 

CVJ 

i 

o 

0. 092 + . 024 

0. 108 + . 014 

0. 114 +. 011 

0. 108 + . 023 

0.125 

2. 0 - 3. 7 

0. 058 + . 015 

0. 058 + . 008 

0. 057 + . 005 

0. 052 + . Oil 

0. 050 

3. 7 - 6. 0 

0. 041 +.011 

0. 037 + . 005 

0. 032 + . 003 

0. 023 + . 005 

0. 034 



Rate 

(counts /cm^- sec) 



0. 6 - 6. 0 

0. 334 + . 097 

0. 358 + . 060 

0. 365 + . 049 

0. 347 + . 088 

0. 364 



Brecher- Morrison Background Subtracted 
Rate (counts /cm^- sec- Me V) 


Energy Loss 
Range (MeV) 

Z = 50 
a rel = - 21 

Z = 70 

a = . 067 

rel 

Z = 100 

a =.15 

rel 

Z = 150 

a = . 36 

rel 

Obs. 

Excess 

0. 6 - 1.0 

0. 110 + . 021 

0. 140 + .011 

0. 165 + . 022 

0. 191 + . 059 

0. 141 

1 . 0 - 2. 0 

0. 081 + . 015 

0. 091 + . 007 

0. 093 + . 012 

0. 084 + . 026 

0. 103 

2. 0 - 3. 7 

0. 051 + . 010 

0. 049 + . 004 

0. 046 + . 006 

0. 040 + . 012 

0. 044 

3. 7 - 6. 0 

0. 036 + . 007 

0. 031 + . 002 

0. 026 + .003 

0. 018 + . 006 

0. 031 



Rate (counts /cm^- sec) 



0. 6 - 6. 0 

0.295 + . 056 

0. 302 + . 023 

0. 297 + . 040 

0.269 + . 083 

0. 308 
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LIST OF FIGURES 


1. The spectrum measured in interplanetary space with the anti- 
coincidence sheild turned on and off. The >25keV rate has been 
corrected to differential flux by letting AE = 225 keV. The differ- 
ence is compared with that predicted from cosmic- ray energy losses 
in the crystal edges. The anti-on rates are interpreted as due to 
cosmic -y-rays. 

2. Data obtained from various spacecraft using omniu. /ectional Nal 
scintillation counters on the total cosmic y-ray flux. Measured 
energy-loss spectra above 200 keV require unfolding of detector 
response for comparison with theoretical spectra; such unfolding 
necessarily involves assumptions about the photon fluxes at all 
energies above 200 keV. 

3. Theoretical model spectra which fit the ERS-18 measured data on 
cosmic y-ra y fluxes. The two models for the isotropic background 
flux are adjusted to the data in the 30 - 300 keV range. The Stecker 
model spectra are adjusted to provide the best fits to the ERS-18 
excess above the two background models. The total photon spectrum 
obtained by summing the background and Stecker models is also 
shown. 
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